One-neutron removal strength of 7 He into 6 He using the complex scaling method 
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We study the one-neutron removal strength of the 7 He ground state, which provides us with 
the 6 He-n component in 7 He. The He isotopes are described on the basis of the 4 He+Xn cluster 
model (X — 1,2,3). The complex scaling method is applied to describe not only the Gamow 
resonances but also the nonresonant continuum states of valence neutrons, with the correct boundary 
condition of particle decays. The one-neutron removal strength of 7 He into the unbound states of 
6 He is calculated using the complex-scaled Green's function, in which a complex-scaled complete 
set of 4 He+n+n states is adopted. Using this framework, we investigate resonant and nonresonant 
contributions of the strength, which individually produce specific structures in the distributions. 
In addition, we propose a method to obtain the real-value strength using the complex values of 
spectroscopic factors of Gamow states. As a result, the 6 He(2 + ) resonance is found to give the 
largest contribution. 

PACS numbers: 21.60.Gx, 21.10.Pc, 27.20. +n 
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I. INTRODUCTION 

The development of experiments using radioactive 
beam has provided us with much information on unstable 
nuclei far from the stability. In particular, the light nu- 
clei near the drip-line exhibit new phenomena of nuclear 
structures, such as the neutron halo structure in 6 He and 
"Li [1, 2]. 

Recently, many experiments [3-11] on 7 He have been 
reported and have confirmed that its ground state is as- 
signed to be the 3/2~ resonant state at 0.3~0.5 MeV 
above the 6 He+n threshold energy. For the excited 
states, 7 He can decay not only to two-body 6 He+n chan- 
nels, but also to many-body channels of 5 He+2?i and 
4 He+3n, because 6 He is a Borromean nucleus and breaks 
up easily into Hc+n+n. This multiparticle decay condi- 
tion makes it difficult to settle the excited states of 7 He. 
In fact, contradictions in the observed energy levels and 
their spins still remain. The excited state at E x ~ 3 
MeV has been assigned as 5/2~ in several experiments 
[3, 4, 6, 8, 9]. The existence of l/2~ is also expected 
[5-9] , but the energies and decay width are not yet fixed. 
This l/2~ state with a possibility of the LS partner of 
the ground state is of interest because the LS splitting 
energy in 7 He gives information on the LS interaction in 
drip- line nuclei. Some experiments [5, 8, 9] report 1/2 - 
at around 1 MeV in the low excitation energy. However, 
other observations [6, 7] suggest a little higher excitation 
energy. The spectroscopic factor associated with the 6 He 
halo state has also been reported [10]. 

On the theoretical side, ab initio calculations of the 
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no-core shell model [12] and the Green's function Monte 
Cairo [13] have been performed, and the calculated en- 
ergies of the ground and 5/2~ states show a good corre- 
spondence with the experiments. The l/2 _ state is pre- 
dicted at around 3 MeV, although the results depend on 
the choice of the three-nucleon forces [13]. Those calcu- 
lations are based on the bound state approximation and 
the continuum effect from many-body open channels is 
not included, while the observed 7 He states are unbound. 

Several methods have been proposed to treat the con- 
tinuum effects explicitly, such as the continuum shell 
model [14] and the Gamow shell model [15-17]. It is, 
however, difficult to satisfy the multiparticle decay con- 
ditions correctly for all open channels. The obtained en- 
ergy spectra and decay widths of He isotopes depend on 
the treatment of open channels. For the spectroscopy of 
7 He, it is necessary to describe the four-body resonances 
of the 4 He+3n system in the theoretical model. Further- 
more, it is important to reproduce the threshold ener- 
gies of the particle decays into subsystems. Emphasizing 
these theoretical conditions, in our previous work [18], 
we investigated the He spectroscopy with the appropri- 
ate treatment of the decay properties concerned with the 
subsystems of 5 ' 6 He. We employed the cluster orbital 
shell model (COSM) of 4 Hc+ln [19-21], in which the 
open channel effects of the 6 He+n, 5 He+2n and 4 He+3n 
decays are taken into account explicitly. This means 
that our analysis satisfies the simultaneous descriptions 
of 4 ' 5,6 He. We described the many-body resonances by 
using the complex scaling method (CSM) [22-25], under 
the correct boundary conditions for all decay channels. In 
CSM, the energies and decay widths of many-body res- 
onances are directly obtained by diagonalization of the 
complex-scaled Hamiltonian with L 2 basis functions. 

We found five resonances of 7 He and investigated their 
properties, such as energies, decay widths, configurations 
and spectroscopic factors (S factors) of the 6 He-n com- 



ponents. In particular, S factors are important in under- 
standing the coupling between 6 He and a last neutron in 
7 He. In the analysis of the S factors of the Gamow states 
of 6 ' 7 He, the S factors are often described in complex val- 
ues, because Gamow states have complex eigenvalues and 
their matrix elements are also complex values [25, 26]. 
These complex S factors lead to the problem of physical 
interpretation with respect to the observation, which was 
discussed in a previous study [18]. To avoid this prob- 
lem, in this article, we develop the method of considering 
the contributions of nonresonant continuum states for £ 
factors in addition to those of Gamow states. This is 
performed by employing the complete set consisting of 
Gamow resonant and nonresonant continuum states in 
the continuum energy region. Using this complete set of 
unbound states of 6 He, we show the one-neutron removal 
strength of the 7 He ground state into the 6 He unbound 
states. This strength function reflects the 6 He-n compo- 
nents of 7 Hc and is observable. 

We explain the relation between the one-neutron re- 
moval strength and the S factors. In Fig. 1, we de- 
pict the schematic illustration of the one-nucleon removal 
strength S(E) from a mass A system to an A — 1 one. 
The strength function S(E) involves the information of 
bound, resonant, and nonresonant continuum states of 
the A — 1 final states. The resonant parts of the strength 
correspond to the S factors of resonances. Beyond the 
threshold energy E t h , the unbound states start to appear, 
and the strength becomes continuous. It is, however, dif- 
ficult to describe the unbound states beyond the two- 
body case, like 4 Hc+n+n of 6 Hc. In the previous study 
[18], only the resonant component (S factors) of the fi- 
nal states is considered as a part of S(E), and the non- 
resonant component is missing. Because of this partial 
treatment of the unbound states, the relation between 
the complex S factors and the real observable S(E) is 
obscure. Hence, in this study, we include the contribu- 
tions not only of the resonances, but also of the nonreso- 
nant continuum states in the distribution, S(E) in Fig. 1. 
Some resonances of the A — 1 system can make structures 
in S(E) like those shown in Fig. 1. It is interesting to 
clarify the origins of the structures seen in the strength. 

In this article, we calculate the one-neutron removal 
strength of 7 He by using the complex-scaled 4 Hc+n+ra 
complete set of the 6 He final states [27], which consists 
of not only the 6 He resonances, but also the contin- 
uum states of 5 He+n and 4 He+n+n. We express the 
6 He complete set by using the complex-scaled solutions 
of the 4 He+n+ra model. It has been shown that CSM 
is a powerful method to investigate the resonances and 
the nonresonant continuum states. So far, using CSM, 
the Coulomb breakup strengths of halo nuclei have been 
successfully analyzed [25, 27-31]. In calculations of the 
strength S(E) in this study, we adopt the complex-scaled 
Green's function described by using the complex-scaled 
eigenstates of 6 He. Similar methods using complex- 
scaled Green's function have been used successfully in 
solving the nuclear reaction problems [32-34]. 
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FIG. 1: Schematic illustration of the one-nucleon removal 
strength S(E), from the mass A system to the A — 1 one. 
The excitation energy E is for the A — 1 system. E t h is the 
lowest threshold energy of particle emission for the A — 1 sys- 
tem. 

In this study, we evaluate the one-neutron removal 
strength of 7 He into 6 He as an observable. This means 
that we develop a method to obtain the observable from 
the complex S factors of Gamow states. Furthermore, we 
investigate the structures in the strength by decomposing 
the 6 He unbound states into 6 He resonances, 5 He+n and 
4 He+n+n continuum components [27]. As was shown 
in the Coulomb breakups of halo nuclei, the contribu- 
tions of each component can be unambiguously classified 
by using CSM. This is a prominent point of the present 
method. 

In Sec. II, we explain the complex-scaled COSM wave 
function for He isotopes, and the method of calculating 
the neutron removal strength using CSM by introduc- 
ing the complex-scaled Green's function. In Sec. Ill, we 
discuss the 7 He structures, the £ factors of the 6 He-n 
components, and also the one-neutron removal strength 
into 6 He. A summary is given in Sec. IV. 

II. COMPLEX-SCALED 4 HE+XiV COSM FOR 
HE ISOTOPES 

A. COSM for the 4 He+In systems 

We explain the COSM of the 4 He+A"n systems, in 
which X = 1 for 5 Hc, X = 2 for 6 Hc and X = 3 
for 7 Hc. The Hamiltonian is the same as that used in 
Refs. [18, 27]; 

X + l X X 

i— 1 i— 1 i<j 

where ti and Tq are the kinetic energies of each parti- 
cle (Xn and 4 He) and of the center of mass of the total 
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FIG. 2: Sets of the spatial coordinates in the COSM for the 
4 He+Xn system. 



system, respectively. The interactions V an and V nn are 
given by the modified KKNN potential [35] and the Min- 
nesota potential [36], respectively. They reproduce the 
low-energy scattering data of the 4 He-n and the n-n sys- 
tems, respectively. 

For the wave function, 4 He is treated as the (Os) 4 con- 
figuration of a harmonic oscillator wave function, whose 
length parameter b c is 1.4 fm to fit the charge radius 
of 4 He. The motion of valence neutrons around 4 He is 
solved variationally using the few-body technique. We 
expand the relative wave functions of the He+Xn sys- 
tem using the COSM basis states [19-21]. In COSM, the 
total wave function ^ of the 4 He+Xn system is given by 
the superposition of the configuration ^> c as 

*( 4 He + Xn) = ^C c * c ( 4 Hc + Xn), (2) 

c 
X 

M/ c ( 4 Hc + Xn) = JJo+jO), (3) 



i=i 

where 4 He corresponds to a vacuum |0). The creation op- 
erator ajj. is for the valence neutron above 4 He, with the 
quantum number Qj in a jj coupling scheme. The index 
i = 1, • • • , X is for X valence neutrons. The set of on is 
included in the index c. The variational coefficient is ex- 
pressed by C c with respect to <I/ C . We take a summation 
over the available configurations. The coordinate repre- 
sentation of the single-particle state corresponding to a^ a . 
is given as ip ai with the relative coordinate r, between 
4 He and a valence neutron shown in Fig. 2. Considering 
the angular momentum coupling, the total wave function 
\E' J with the spin J is expressed as 



$ J ( 4 Hc + ln) 
* c J ( 4 Hc + ln) 



]Tc/* c J ( 4 Hc + An), 



^{[$( 4 Hc), Xc J (Xn)] J } 
X J c {n) = 
xf(2n) = A{[ip ai ,ip a2 ]j}, 
xt(3n) = A{[[^ ai ,^ a2 } jl2 ,^ a3 ]j}. 



(4) 

(5) 
(6) 
(7) 
(8) 



Here, as shown in Fig. 2, \t (Xn) expresses the wave 
functions for valence neutrons. The spin ji2 is a coupled 



angular momentum of the first and second valence neu- 
trons. The antisymmctrizcrs between valence neutrons 
and between a valence neutron and neutrons in 4 He are 
expressed as A and A', respectively. The effect of A' is 
treated in the orthogonality condition model [21, 25, 27], 
in which ip a is imposed to be orthogonal to the Os state 
occupied by 4 He. The radial part of tp a is expanded with 
a finite number of Gaussian basis functions as 

i> a = J2 Ca > k <^( r > 6 ^fc)> ( 9 ) 
fc=i 

c/> k a (r,b a , k ) = Mr^e-^^^YeJ^^i^-m 

Here, the index k is for the Gaussian basis with the length 
parameter fe Q .fc. A basis number for the state a and the 
normalization factor for the basis are given by N a and 
A", respectively. Two expansion coefficients {C/} and 
{C a ,k} arc determined variationally with respect to the 
total wave function ^> J by diagonalization of the Hamil- 
tonian matrix elements. The length parameters b ai k are 
chosen as geometric progression [25]. Due to this expan- 
sion of the radial components using a finite number of 
basis states, all the energy eigenvalues are discretized for 
bound, resonant and continuum states. We use at most 
17 Gaussian basis functions with a maximum range of 40 
fm. 

Here, we discuss the coupling between 4 He and va- 
lence neutrons. This is related to the boundary condi- 
tion of the neutron emission, which is important, in par- 
ticular, when we deal with the weakly binding, resonant 
and continuum states [18, 31]. We consider the coupling 
between 7 He and the 6 He+n configurations. Asymptoti- 
cally, when the last neutron is located far away from 6 He 
(r 3 — > oo for 4 He+3n in Fig. 2), any coupling between 
6 He and a last neutron disappears, and 6 He becomes its 
isolated eigenstate of the Hamiltonian in Eq. (1). 



* J ( 7 He) 



C c J ^'{[$( 4 He), Xc J (3n)] J } (11) 



*^'( 6 He) 



[<( 6 He),V> Q3 ] J , (12) 
^^>'{[$( 4 Hc), Xc ^(2n)] J '} 1 3) 



where the spin J and J' are for 7 He and 6 He, respec- 
tively, and the index v indicates the eigenstate of 6 Hc. 
The mixing coefficient and the two-neutron wave 

function \i v{^ n ) m Eq. (13) are those of the 6 He eigen- 
states. Hence, the three- neutron wave function %^(3n) 
in Eq. (11) satisfies the following asymptotic forms 



c \ c 



n) Ua,. (14) 



This relation implies that the three-neutron wave func- 
tion of 7 He is asymptotically decomposed into two- 
neutron wave function of 6 He and a last neutron. Equa- 
tions (11)-(14) describe the boundary conditions of the 
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COSM wave functions for He isotopes. In contrast, when 
a last neutron comes close to 6 He, this last neutron dy- 
namically couples to the 6 He eigenstates ^ J v . This cou- 
pling depends on the relative distance between 6 He and a 
last neutron, and changes the 6 Hc+n configurations from 
the eigenstates of 6 He to the eigenstates of 7 He. In the 
COSM, the structure change of 6 He inside 7 He is deter- 
mined variationally to optimize the 7 He eigenstate. This 
discussion of asymptotic condition can be also applied 
for the configurations of 5 He+2n and 4 He+3n. Hence, in 
the COSM, we can treat the neutron emissions with the 
correct boundary conditions. 

For the single-particle states a = £j (j = I <g> 5), we 
take angular momenta I < 2 to keep the accuracy of the 
converged energy within 0.3 MeV. In the calculation of 
7 He, we adjust the calculated energies of 6 He by taking 
the 178.8 MeV of the repulsive strength of the Minnesota 
force [36] and the three-cluster interaction V ann for the 
4 He-n-n system [27]. The former adjustment of the NN 
interaction can be understood from the pairing correla- 
tion between valence neutrons with higher angular mo- 
menta I > 2 [35]. The latter is considered to come from 
dominantly the tensor correlation in 4 He. Recently, we 
showed that the binding energy and the excited states of 
6 He can be well explained without the three-cluster in- 
teraction by taking into account the tensor correlation 
of 4 He explicitly [37, 38]. Following previous studies 
[18, 27], we use the three-cluster potential: 

yann = e -tf+r])/bl with V3 = _25 MeV.(15) 

i<j 

Adding this potential to the Hamiltonian in Eq. (1), our 
results agree with the observed energies of 6 Hc for + 
and 2 + states. Hence, the present model reproduces the 
observed properties of 5 ' 6 He, simultaneously [18, 39], as 
shown in the next section (Fig. 4), namely, the threshold 
energies of the particle emissions from 7 He. 



B. Complex scaling method (CSM) 

We explain the CSM, which describes resonances and 
nonresonant continuum states. Hereafter, we refer to the 
nonresonant continuum states as simply the continuum 
states. In the CSM, we transform the relative coordinates 
of the 4 Ke+Xn model shown in Fig. 2, by the operator 
Ug as 

U e : r t r,e ! " for i = 1, • • • , X , (16) 

where 9 is a scaling angle. The Hamiltonian in Eq. (1) 
is transformed into the complex-scaled Hamiltonian 
Hg = UgHUg 1 , and the corresponding complex-scaled 
Schrodinger equation is given as 

H e ^i = E^ J , (17) 

= e (3/2)i0-X^7 ( | r . e « }) (lg) 



The eigenstates tyg are obtained by solving the eigen- 
value problem of Hg in Eq. (17). In CSM, we obtain all 
the energy eigenvalues E of bound and unbound states on 
a complex energy plane, governed by the ABC theorem 
[22] . In this theorem, it is proven that the boundary con- 
dition of Gamow resonances is transformed to the damp- 
ing behavior at the asymptotic region. This condition 
enables us to use the same theoretical method to obtain 
the many-body resonances as that used for the bound 
states. For a finite value of 9, every Riemann branch cut 
is commonly rotated down by 29. For 7 He, the contin- 
uum states of 6 He+n 5 He+2n and 4 He+3n channels are 
obtained on the branch cuts rotated with the 29 depen- 
dence [18]. On the contrary, bound states and resonances 
are discrete and obtainable independently of 9 (see Fig. 
3 for the 6 Hc(2 + ) case). Hence, these discrete states are 
located separately from the many-body continuum spec- 
tra on the complex energy plane. We can identify the 
resonance poles of complex eigenvalues: E = E r — iT/2, 
where E r and T are the resonance energies measured from 
the lowest threshold and the decay widths, respectively. 

In the wave function, the 9 dependence is included in 
the variational coefficients in Eqs. (4) and (9) as {C/' e } 
and {C^ fc }, respectively. We take the value of 9 as 30° 
in this calculation. 

In the study, we calculate the one-neutron removal 
strength of He into 6 He, in which we need a complete set 
of 6 Hc including bound, resonant, and continuum states. 
We express this 6 He complete set using the complex- 
scaled eigenstates tyg obtained in the 4 He+n+n model. 
Hereafter, we denote the 6 He wave function with the state 
v' as simply $„/, and the 7 He wave function ^^C^He) as 

We briefly explain the extended completeness relation 
(ECR) of 6 He using CSM [27, 40, 41]. When we take a 
large 9 like in Fig. 3, three-body unbound states of 6 He is 
decomposed into three categories of discrete three-body 
resonances, three-body continuum states of 4 Hc+n+n, 
and two-body continuum states of 5 He(3/2~, l/2~)+n. 
Here, the 5 He+n two-body continuum states are ob- 
tained on the branch cuts, whose origins are resonance 
positions of 5 Hc(3/2~, l/2~), as shown in Fig. 3. Using 
all the unbound states of 6 He, we introduce the extended 
three-body completeness relation (ECR) of the complex- 
scaled Hamiltonian Hg of 6 He as 

1 = tjom 

"V 

= {three-body bound state of 6 He} 

+ {three-body resonance of 6 He } 

+ {three-body continuum states of 4 He+n+n} 

+ {two-body continuum states of 5 Hc+n} , (19) 

where {$^,$^1 are the complex-scaled 6 He wave func- 
tions and form a set of biorthogonal bases. This relation 
is an extension of the two-body ECR [25, 40]. Because 
the detailed definition of the biorthogonal bases is writ- 
ten in the previous works [27, 40], we briefly explain it 
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here. When the wave number k v of <!>„ is for discrete 
bound and resonance states, the adjoint wave number 
k v of is defined as k v = — k*, which leads to the re- 
lation <f>„ = ($„)* [40-42]. For continuum states, the 
same relation of the biorthogonal states of resonances 
is adopted, because we use a discretized representation. 
In the 4 He+n+ra model, the 4 He-2n channel is included 
in the three-body continuum components of Hc+ra+n, 
because 2n does not have any bound states or physical 
resonances. 



and of the creation and annihilation operators of valence 
neutrons, respectively. It is found that the summed value 
of S factors for the 6 He states becomes the number of 
valence neutrons in 7 He. This property of the S factors 
is also established when the complex scaling is operated. 

In the numerical calculation, we express the radial part 
of the operator a a in Eq. (21) using a complete set ex- 
panded by 40 Gaussian basis functions with the maxi- 
mum range of 80 fm for each orbit. This treatment is 
sufficient to converge the results. 



C. Spectroscopic factor of 7 He 



D. One-neutron removal strength of 7 He 



We explain the S factors of the 6 He-n components for 
7 He. This S factor is used in the calculation of the one- 
neutron removal strength of 7 He. As was explained in a 
previous study [18], because Gamow states generally give 
complex matrix elements, S factors of Gamow states are 
not necessarily positive definite and are defined by the 
squared matrix elements using the biorthogonal property 
[16, 18, 27, 41] as 



a 

l 



j>y, a - 2J+1 



(*t'\K\\*t 



J\2 



(20) 
(21) 



where the annihilation operator a a is for valence neutron 
with the state a. The spin J and J' are for 7 He and 6 Hc, 
respectively. The index v (z/) indicates the eigenstate of 
7 He ( 6 He). In this expression, the S factors Sjf v , arc 
allowed to be complex values. In general, an imaginary 
part of 5* factor often becomes large relative to the real 
part for resonances having large decay widths. 

The sum rule value of S factors, which includes Gamow 
state contributions of the final states, can be considered 
[18, 25]. When we count all the obtained complex £ 
factors not only of the Gamow states but also of the con- 
tinuum states in the final states, the summed value of S 
factors satisfies the associated particle number, which is 
a real value. In such a case, the imaginary part of the 
summed S factors becomes zero, similar to the transition 
strength [25, 27, 29]. For 7 He into the 6 He-n decomposi- 
tion, the summed value of the S factors Sjf v , in Eq. (21) 
by taking all the 6 He states, is given as 

E s ry = E ^l M \< m ^A*l M ) 

J' ,v' a,m 

= 3 , (22) 
where we use the completeness relation of 6 He as 



1 = E t)* J J M ')(H M '\. 

J'M' V 



(23) 



Here the labels, M (M') and m are the z components of 
the angular momenta of the wave functions of 7 He ( 6 Hc) 



We explain the one-neutron removal strength of 7 Hc 
into 6 Hc. This is a function of the real energy of 6 Hc, 
E. We first introduce the Green's function G(E, t], rj') of 
6 He, which is used in the derivation of the strength. The 
procedure is the same as that used in the case of the elec- 
tric transitions [27, 29]. The coordinates, rj and rj' , rep- 
resent the set of ?"j (i = 1, ■ • ■ , X) in Fig. 2. Here, we in- 
troduce the complex-scaled Green's function Q 9 {E 1 rj, rj') 
of 6 He as 



G 9 (E,r,,r,') 



E-H 
1 



(24) 



-i 



\E-Hg 

few 

E-Et 



(25) 



In the derivation from Eq. (24) to Eq. (25), we insert the 
ECR of 6 He given in Eq. (19). The 6 He energy, E e v , corre- 
sponds to the eigen wave function <fr e u . The 9 dependence 
of E e v appears only in the continuum spectra. 

Next, the strength function S(E) for the annihilation 
operator a a is defined using the Green's function in a 
usual case without CSM as 

S(E) = ^S a (E), (26) 

S «( E ) = i<*ok|$„)($„|a a |* > 6{E-E V ) 
■v 



x g(E, n^aMv') 



(27) 



For simplicity, we omit the labels of the angular momenta 
and their z components of the wave functions and of the 
operators. The wave function is the ground state of 
7 He. We also consider the sum rule value of the strength 
S(E) in Eq. (26), which is defined by the integration of 
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S{E) over the real energy E of 6 He. Using the complete- 
ness relation of the final states of 6 He, the sum rule value 
is given as 

dES(E) = £^(§ |4l^)(^KI*o) 

a "V 

= X)<*o|4oa|*0> 



= 3. 



(28) 



Thus, it is also confirmed that the energy integrated value 
of S(E) satisfies the associated particle number of 7 Hc, 
similar to the case shown in Eq. (22). 

To calculate the strength function S a (E) in Eq. (27), 
we operate the complex scaling on S a {E), and use the 
complex-scaled Green's function of Eq. (25) as 



S a {E) = — Im 

7T 



x g^E^^ayiW) 
= £s a ,„(E), (29) 

<*ol(4) fl |*v><^l<*al*0 



S au {E) = Im 

7T 



E-Ei 



(30) 



In Eq. (30), the strength function is calculated using 

the one-neutron removal matrix elements (^la^l^o)- It 
is noted that the function S atV {E) is independent of 9 
[18, 27, 40, 43]. This is because any matrix elements 
are obtained independently of 6 in the complex scaling 
method, and also because the state v of 6 He is uniquely 
classified according to ECR defined in Eq. (19). As a 
result, the decomposed strength S a ,v(E) is uniquely ob- 
tained. Thus, the one-neutron removal strength S a (E) 
is calculated as a function of the real energy E of 6 He. 
When we discuss the structures appearing in S a (E) , it is 
useful to decompose S a (E) into each component S atV {E) 
by using the complete set of the final state v of 6 He. We 
can categorize v of 6 He using the ECR in Eq. (19). Be- 
cause of this decomposition of unbound states, we can 
unambiguously investigate how much each resonant and 
continuum state of 6 He exhausts the strength. This is a 
prominent point of the present method. 

Here, we discuss the properties of the decomposed 
strength function S a , v {E) in Eq. (30). The original 
strength function S a (E) in Eq. (29) corresponds to 
an observable being positive definite for every energy. 
However, the decomposed strength function S a ^{E) is 
not necessarily positive definite at all energies, because 
S a ,v{E) cannot be directly observed, similar to the res- 
onant poles. This means that S ayV (E) can sometimes 
have negative values. This property of the decomposed 
strength has been generally discussed in the electric tran- 
sition strength [27, 29]. 

The stability of the calculated matrix elements of res- 
onant and continuum states using the CSM has been 



shown in many works [25, 27, 40]. For continuum states, 
we adopt the discretized representation using the L 2 in- 
tegrable basis functions. This discretization has been 
checked to reproduce the genuine continuum states by 
using the CSM [32, 33, 43]. 



III. RESULTS 
A. Energy spectra of He isotopes 

We first discuss the energy spectra of 5_7 He obtained 
in the present model. The example of the eigenvalue 
distribution of 6 He(2 + ) obtained in the CSM is shown 
in Fig. 3. In this figure, we can easily identify the loca- 
tions of the three-body resonances of 6 He (2+ and 2^~) 
and, further, of the two kinds of continuum states of 
5 Hc(3/2~,l/2~)+n and 4 He+n+n. In Fig. 4, we summa- 
rize the energy spectra of 5_7 He. A detailed discussion 
of the obtained energy levels is given in previous works 
[18, 25]; and hence we briefly explain it here. 

In Tables I and II, we list the energies and decay widths 
of 6 He and 7 He, respectively, measured from the 4 He+Xn 
threshold energy. For 7 He, the ground state energy E r is 
0.184 MeV measured from the 6 He+n threshold, which 
is slightly overbound with respect to the experiments 
{E r = 0.44(2) MeV [3] and 0.36(5) MeV [8]). Due to 
this overbinding, the calculated decay width is smaller 
than the experiments of T ~ 0.16 MeV [3, 8]. When 
we fit the observed energy of E r = 0.44 MeV, namely, 
—0.54 MeV from the 4 He+3n threshold, by reducing the 
strength of V ann , the decay width becomes 0.14 MeV 
and nicely agrees with the experiments, as is shown in 
Table II. 

We found the 5/2~ state, whose position agrees with 
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0.5 1 1.5 2 2.5 3 3.5 4 
Re(E) [MeV] 



FIG. 3: (Color online) Energy eigenvalues of the 6 He(2 + ) 
states in the complex energy plane measured from the 
4 He+n+n threshold. Three schematic lines with dots cor- 
respond to the 4 He+n+n and 5 He(3/2~ ,l/2~ )+n continuum 
states, in order from left to right, respectively. 
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FIG. 4: Energy eigenvalues of the obtained 5 ' 6 ' 7 He resonances 
measured from the 4 He+Xn threshold. 



TABLE I: Energy eigenvalues of the He resonances measured 
from the 4 He+2n threshold. The values with parentheses are 
the experimental ones [39] . Dominant configurations are also 
listed. 





Energy (MeV) 


Width (MeV) 


Config. 




-0.974 (-0.975) 




O3/2) 2 


o 2 + 


3.90 


8.77 


(Pi/2) 2 


2f 


0.840 (0.822(25)) 


0.107 (0.113(20)) 


(P3/2) 2 


2 2 + 


2.53 


3.87 


(P3/2)(Pl/2) 


1+ 


3.02 


5.97 


(P3/ 2 )(Pl/ 2 ) 



the several experiments [3, 4, 8]. Further, the 3/2,7 state 
is degenerated with the 5/2" state and their decay widths 
do not differ so much. This result suggests the possibility 
of the superposed observation of two states in this energy 
region. We found one broad 1/2 - resonance at a low ex- 
citation energy of E x =1.05 MeV, while the experimental 
uncertainty is still large [5-9, 11]. Recently, there has 
been a report discussing the order of energy levels [11], 
which fully agrees with our results. 

For the neutron configurations, the resonances of 6 Hc 
and 7 He are dominantly described by the p-shell con- 
figurations and the small contributions come from the 
sd-shell [18]. The dominant configurations for each res- 
onance are listed in Tables I for 6 He and Table II for 
7 He, respectively. From the configurations, several ex- 
cited states of 7 He seem to have the 6 He+n type configu- 
rations. To understand the detailed structures of the 7 He 
states, we investigate the S factors of the 6 Hc+n system 
and the one-neutron removal strength of 7 He. 

B. Spectroscopic factors of 7 He 

We calculate the 5* factors of the 6 He-n components of 
the 7 He resonances. Here we choose the and 2+ states 
of 6 He. It is noted that the present S factors correspond 
to the components of 6 He inside the 7 He resonances and 
contain the imaginary part. We improve the operation 
of complex scaling to calculate the matrix elements of 



TABLE II: Energy eigenvalues of the 7 He resonances mea- 
sured from the 4 He+3n threshold. The values with parenthe- 
ses are the ones fitted to the position of the observed reso- 
nance energy of the ground state. Dominant configurations 
are also listed. 





Energy (MeV) 


Width (MeV) Config. 


3/27 


-0.790 (-0.54) 


0.014 (0.14) (p 3/2 ) 3 


3/27 


2.58 


1-95 (P3/2) 2 (Pl/ 2 ) 


3/2 3 - 


4.53 


5.77 (P3/2XP1/2) 2 


1/2- 


0.26 


2.19 (P3/2) 2 (Pl/ 2 ) 


5/2" 


2.46 


1.50 (p 3 /2) 2 (pi/2) 


TABLE III: Spectroscopic factors of the 6 He-n components of 


He. Details are described 


in the text. 




6 He(0+)-n 


6 He(2+)-n 




Present 


CK 


Present CK 


3/27 


0.64-M0.06 


0.59 


1.55-«0.31 1.21 


3/27 


0.005 -M0.01 


0.06 


0.95 + i0.02 1.38 


3/27 


0.003 + iO.0002 




0.02 - i0.004 


1/2- 


1.00 - i0.13 


0.69 


0.10-i0.05 0.60 


5/2" 


0.00 4-i0.00 


0.00 


0.95 + i0.02 1.36 



S factors from the previous ones [18]. In Table III, we 
list the results of S factors of each 7 He resonance. For 
reference, the results of the conventional Cohen-Kurath 
shell model (CK) [6] are shown with real values due to the 
bound state approximation of resonances. With regard 
to the real parts of the calculated S factors, the trend 
seen in our results is roughly similar to that of the CK 
results. In Table IV, the S factors are decomposed into 
the components of two p orbits of the last neutron. It is 
found that every 7 He resonance is dominated by one of 
the p orbits coupled with 6 He. 

In Table III, for the 3/27 state, the 6 He(0^)-n com- 
ponent almost shows a real value with a small imaginary 
part. This real part well corresponds to the recent ob- 
servation of 0.64 ± 0.09 [10]. The 6 He(2+)-n component 
is large, more than twice of that of the 6 Hc(0 ] f )-n case. 
For the 3/2^" state, 6 He(2+) is dominantly mixed, and 
this characteristic is common for the 5/2" state. Hence, 
these two states have a similar structure of the valence 
neutron configurations. For the 3/2.7 state, two compo- 
nents of 6 He(0f , 2f )-n are very small. This state can be 
considered to have some components of 6 He(0^,l + ,2j)-n 
from the configurations shown in Tables I and II. 

For the 1/2" state, the S factor of 6 He(0^)-n(p 1/2 ) 
is almost unity with a small imaginary part, and the 
6 Hc(2^)-n component is small. These results of l/2~ 
indicate that the 6 He(0^)-ra component is dominant in 
7 He(l/2") and also suggest the weak coupling nature of 
the P1/2 orbital neutron around 6 He, which retains two- 
neutron halo structure. Therefore, the l/2~ resonance 
could be mainly considered as a single particle resonance 
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TABLE IV: Components of two p orbits of the last neutron 
in S factors of 6 He-n(p-wave) of 7 He. Upper and lower tables 
correspond to 6 He(0^) and 6 He(2^), respectively. 





6 He(0+)-p 3 /2 


6 He(0+)-p 1/2 


3/2r 


0.64 + i0.06 


0.00 + i0.00 


3/2 2 


0.005 -M0.01 


0.00 + iO.OO 


3/2 3 - 


0.003 + iO.0002 


0.00 + iO.OO 


1 IT~ 


n no -+- j'o no 


1 DO — iO 1 3 


5/2 _ 


O.OO + iO.00 


0.00 + iO.OO 




6 He(2+)-p 3/2 


6 He(2+)-p 1/2 


3/2" 


1.54-i0.31 


0.005 - t0.002 


3/2 2 " 


0.001 + iO.001 


0.95 + i0.02 


3/2 3 " 


0.007 - i0.015 


0.02 + i0.01 


1/2" 


0.10-i0.05 


O.OO + iO.OO 


5/2" 


0.10 + i0.02 


0.85 - iO.001 



of the pi/2 neutron surrounding a halo state of 6 He. 

We obtained the interesting information of the struc- 
tures of the 7 He resonances via S factors. We should be, 
however, careful to derive the conclusion of physical in- 
terpretation of the complex S factors of Gamow states, as 
was mentioned in the previous studies [18, 27, 29]. In this 
analysis, both of the initial ( 7 He) and final ( 6 He) states 
are Gamow states. In the next subsection, we discuss a 
one-neutron removal strength from the 7 He ground state 
assuming a bound state wave function to the continuum 
energy states of 6 He. The continuum energy states of 
6 He are described using the complex-scaled complete set 
of 6 Hc. It is shown that we can see contributions not only 
of resonances but also of nonresonant continuum states 
of 6 Hc to the strength function. 



C. One-neutron removal strength of 7 He 

In Sec. IIIB, the obtained S factors are complex val- 
ues, which are caused from the Gamow states of 7 He 
and 6 Hc(2 ] f ). For the final states of 6 He, we consider 
only resonances, not the continuum states of 5 He+n and 
4 He+n+n. The decay contributions from 7 He to these 
continuum states arc not clarified. In this study, we 
take up this problem and show the strength function 
of the 7 He ground state into 6 He using Eqs. (29) and 
(30). For 6 He, we prepare the three-body complete set 
of the 4 He+n+n model, namely, ECR in Eq. (19) by us- 
ing the complex-scaled wave functions <fr e u of 6 He. In this 
method, we take care not only of the 6 He resonances, but 
also of the continuum states of 5 He+n and 4 He+ra+n. 
Here, we use the bound state approximation for the ini- 
tial 7 He ground state, because this state has a very small 
decay width, as shown in Table II. 

In Fig. 5, we show the one-neutron removal strength 
S{E) of 7 He(3/2~) into 6 Hc with spin J, where the en- 



1 


1 

1 6 He(2|) 


(a) ; 


- 6 He(0{) 
- / 

1 i / 


V , , 





-1 fl 2 3 4 5 
2\ E [MeV] 




1 2 3 4 5 
E [MeV] 

FIG. 5: One-neutron removal strength of 7 He into 6 He(J 7r ,_E) 
states measured from the 4 He+n+n threshold energy. The 
vertical arrow in (a) indicates the resonance energy of 
6 He(2 ] f ). The strength into 6 He(0^~) is shown by a histogram 
at the position of the corresponding energy. 

ergy E of 6 He is measured from the 4 He+n+n threshold. 
We show the results of only plus parity states of 6 He, 
while the negative parity strengths are found to give a 
negligible contribution in the order of 10~ 3 of S(E). In 
Fig. 5, the obtained strengths remain positive definite. 
The dominant component comes from the 2 + state in 
Fig. 5 (a), whose strength has a peak at the resonance 
energy, 0.84 MeV, of 6 Hc(2 ] f ). Below the 4 He+n+n 
threshold, 6 Hc(0 ] f ) gives a contribution of 0.61 of S(E), 
which is close to the real part of 0.64 in Table III. The 
small difference comes from the bound state approxima- 
tion of the He ground state. For the continuum energy 
region above the 4 He+n+n threshold, the + , 1 + and 
3 + states give small contributions in the strengths, as 
shown in Fig. 5 (b). It is also found that three kinds 
of strengths commonly show the enhancement at around 
1 MeV. From these results, it is concluded that the onc- 
neutron removal strength of 7 Hc is dominantly exhausted 
by the 6 He(2^) resonance above the 4 He+n+n threshold 
energy. 

We discuss the detailed structures seen in the four 
strengths in Fig. 5. The origin of the 2 + peak is the 
6 Hc(2 ] f ) resonance, whose matrix element (1.55 — z0.31) 
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in Tabic III is so large that the strength function S(E) 
produces the sharp peak due to the small decay width of 
0.11 MeV of 6 He(2+ ). The remaining continuum strength 
of the 5 He+n and 4 Hc+n+n components are very small 
and masked by the sharp peak from resonance. For the 
1 + and 3 + states, we see one peak and a shoulder-like 
structure. For the + state, only one peak is seen. It is 
interesting to examine the origins of these structures in 
each strength. We do this it by decomposing the strength 
functions using the ECR of 6 He in Eq. (30). 

In Fig. 6, the 2+ strength distribution is decomposed 
into the three kinds of the components of 6 He(2^~) res- 
onance, 5 He(3/2~)+ra and 4 Hc+n+n continuum states. 
It is clearly and explicitly shown that the sharp peak 
is exhausted by the 6 He(2^) resonance. The residual 
2 + continuum strengths are found to be small. Among 
them, the 5 He(3/2~)+n two-body continuum component 
makes a peak at around 0.75 MeV. This energy coincides 
with the position of the 5 He+n threshold energy (0.74 
MeV), and then the peak reflects the threshold effect of 
the 5 Hc+ra open channel. The 4 He+n-|-n contribution is 
small and does not produce a definite structure in the 
distribution. 

In Fig. 7, the strengths of 6 He(0 + ,l + ,3 + ) are decom- 
posed into two kinds of continuum components. In ev- 
ery figure, it is found that the low-energy peak observed 
commonly at around 1 MeV, comes from the contribu- 
tion of 5 He(3/2 - )+n. This peak can be considered as a 
threshold effect of the two-body 5 He(3/2~)+n channel. 
The 4 He+n+n component also commonly makes a mild 
bump at around 1.5 MeV, which makes the shoulder-like 
structure in the total strength for 1 + and 3 + , in Fig. 7 (b) 
and (c), respectively. This Hc+n+n strength is consid- 
ered to correspond to the background contribution and 
does not indicate the existence of any physical states. For 
the + state in Fig. 7 (a), the 4 He+n+n component is 
smaller than the 5 He(3/2~)+n one, and then only one 
peak structure is confirmed in the total strength. As a 



result, it is found that for three spin states, the strengths 
into the 6 He unbound states show some structures, whose 
origins come from the two kinds of different continuum 
components of 5 He+n and 4 He+n+n. 

It is noted that there are other components of broad 
resonances such as 6 He(0^~, 1 + , 2%) and also a binary 
5 He(l/2~)+n one. Their contributions are checked to 
give negligible contributions in the strength. This is be- 
cause the magnitudes of their matrix elements are very 
small, and further, these states have large decay widths 
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FIG. 6: Decomposition of the 6 He(2 + ) component of the 
strength functions. The vertical arrow indicates the threshold 
energy (0.74 MeV) of the 5 He(3/2~)+n channel. 



FIG. 7: Decomposition of the 6 He(0 + ,l + ,3 + ) components of 
the strength functions. The vertical arrows in each panel 
indicate the threshold energy of the 5 He(3/2 - )+n channel. 
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as shown in Table I. For 5 He(l/2~), its decay width is 
obtained as 5.84 MeV. Similar results from broad reso- 
nances were confirmed in the Coulomb breakup strengths 
of halo nuclei [27, 29]. 

It is found in Figs. 6 and 7 that the 5 He(3/2~)+n two- 
body continuum component remains with finite values 
even below the two-body threshold energy. This is be- 
cause the 5 He in the 5 He+n threshold, has a decay width 
of 0.60 MeV. Hence, the 5 He+n component of S(E) can 
have a strength below the threshold energy [27]. 

In this model, the 7 He ground state energy is slightly 
overbound with respect to the experiments. Hence, we 
also calculate the strengths when we fit the observed 
ground state energy of 7 He, as shown in Table II. It 
is confirmed that the trend of the obtained strengths 
does not depend on the 7 He ground state energy, and 
the present conclusion does not change. 

IV. SUMMARY 

We have investigated the structures of 7 He with the 
four-body cluster model. The boundary condition for 
many-body resonances is accurately treated using the 
complex scaling method. In this study, we have devel- 
oped the analysis of 7 He. From the S factors of the 6 He-n 
component of 7 He, the l/2~ resonance has a possibility 
to be the weakly coupled state consisting of a halo state 
of 6 Hc and a pi/ 2 valence neutron. 

We have also proposed the method to obtain the 
real-value strength functions from Gamow states giving 
complex- value S factors. This leads to the calculation of 
the one-neutron removal strength of 7 He into 6 He. We 
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